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Abstract 

Although full wave electromagnetic systems are large and cum- 
bersome to solve, typically only a few parameters, such as input 
impedance, S parameters, and far field pattern, are needed by the 
designer or analyst. A reduced order modeling of these parameters is 
therefore an important consideration in minimizing the CPU require- 
ments. The Asymptotic Waveform E valuation (AWE) method is one 
approach to construct a reduced order model of the input impedance 
or other useful electromagnetic parameters. We demonstrate its appli- 
cation and validity when used in conjunction with the finite element 
method to simulate full wave electromagnetic problems. 


1 Introduction 

The method of Asymptotic Waveform Evaluation (AWE) provides a reduced- 
order model of a linear system and has already been successfully used in VLSI 
and circuit analysis to approximate the transfer function associated with a 
given set of ports/ variables in circuit networks [1, 2, 3]. The basic idea of 
the method is to develop an approximate transfer function of a given linear 
system from a limited set of spectral solutions. Typically, a Pa.de expansion 
of the transfer function is postulated whose coefficients are then determined 
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by matching the Pade representation to tin' available spectral solutions of 
t lie complete system. 

In this paper we investigate the suitability of the AWK method for ap- 
proximating the response of a given parameter in full wave simulations of 
radiation or scattering problems in electromagnetics. Of particular interest 
is the use of AW 1*. for evaluating the input impedance of the antenna over a 
given bandwidth from a knowledge of the full wave solution at a few (even a 
single) frequency points. Also, the method can bo used to fill in a backseat- 
tering pattern with respect to frequency using a few data samples of that 
pattern. Below wo first describe the recasting of the I KM system for ap- 
plication of the AWK. We then proceed to describe the AWE method and 
demonstrate its application, accuracy and efficiency in computing the input 
impedance of a shielded microstrip stub. 


2 Theory of Asymptotic Waveform Evalua- 
tion 

2.1 FEM System 

The application of the finite element method to full wave electromagnetic 
solutions amounts to generating a linear system of equations by extremizing 
the functional [4] 

T =< V xE.o-V xE> —k 2 < EJ3 • E > -ffcb.t. ( 1 ) 

where <. > denotes an inner-product and b.t. is a. possible boundary term 
whose specific form is not required for this discussion. Also, the dyadics A 

and 3 are material related coefficients, k = — = — is the wavenumber and 

A c 

uj is the corresponding operating frequency with c being the speed of light. 
A discretized form of (1) incorporating the appropriate boundary conditions 

is [5] 


(^A 0 + fcAi + k 2 



{.V} = (/} 


( 2 ) 


where A, denote the usual square (sparse) matrices and {/} is a column 
matrix describing the specific excitation. 
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Clearly (2) can lx- solved using direct or iterative methods for a given 
value of the wavenumber. Even though A. is sparse, the solution of the 
system (2) is computationally intensive and must he repeated for each k 
to obtain a frequency response. Also, certain analyses and designs may re- 
quire both temporal and frequency responses placing additional computa- 
tional burdens and a repeated solution of (2) is not an efficient approach in 
generating these responses. An application of AW F. to achieve an approxima- 
tion to these responses is an attractive alternative. Below we formulate AW F. 
in conjunction with the finite' element method (2). for modeling antenna and 
microwave circuit problems. For these problems, t he excit at ion column {/} 
is a linear function of the wavenumber and can therefore he staled as 

{/} = '.{/,} (:t) 

with {/]} being independent of frequency. This observation will be specifi- 
cally used in our subsequent presentation. 

2.2 Asymptotic Waveform Evaluation 

Io describe the basic idea of AWE in conjunct ion with the FEM, we begin 
by first expanding the solution {A'} in a Taylor series about A- 0 as 

{-V} = { A'o} + (k - ko) {A',} + (k - k 0 ) 2 {X 2 } + ... 

+{k — ko) 1 {AJ + O {{k — k 0 ) l+> ] (4) 

where {A' 0 } is the solution of (2) corresponding to the wavenumber k 0 . By 
introducing this expansion into (2) and equating equal powers of k in con- 
junction with (3), after some manipulations, we find that 

{A'o} = AtoAoM/i} 

{.Vi} = Aq 1 [{/i} — A, {A'o} — 2 /c 0 A 2 {A'o}] 

{A 2 } = ~A 0 1 [A] {Aj} + A2 ({Aq} + 2ko {A'] })] (5) 


{A'/} = -Ao 1 [A,{AVi} + A 2 ({AV 2 } + 2A:o{AV 1 })] 

with 


Aq — Aq + k 0 Ai + AtqA 2 


( 6 ) 
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hxprrssions (.’>) are referred to as the system moments whereas |b) is the 
system at 1 Ik* prescribed wavenumber (A‘ 0 ). Although an explicit inversion 
of A " 1 may he needed as indicated in (o). this inversion is used repeatedly 
and can thus he stored out of --con 1 for the' implementation of AW F. Also, 
given that for input impedance computations we are' typically interested in 
the field value at one location of the computational domain, only a single 1 
entry of |A'/(A*)} need he considered, say (the pth entry) Xf'(k). The above* 
moments can then he reduced to scalar form and the expansions (o) hecotnr 
a scalar representation of A 'f(h) about the corresponding solution at h\ } . To 
yield a more* convergent expression, we can instead revert to a Pade expansion 
which is a conventional rational function. 

For transient analysis the Fade expansion can he cast by part ial fract ion 
decomposition [3. 6] into 

■VJ(A •l-.VJ. + tyr^Tfc 

i= 1 

wlicre A ',,0 > s t lif’ limiting value as k tends to infinity. This is a c\th order 
representation and is suitable for time/frequency domain transformation. As 
can be realized, the residues and poles (r, and k 0 + k t ) in (7) correspond 
to those of the original physical system and play important roles in the 
accuracy of the approximation. As can be expected a higher order expansion 
with more zeros and poles can provide an improved approximation. The 
accuracy of AWE relies on the prediction of the dominant residues and poles 
located closest to k 0 in a complex plane. Its key advantage is that for many 
practical electromagnetic problems only a few poles and zeros are needed for 
a sufficiently accurate representation. 

For a hybrid finite element - boundary integral system, the implementa- 
tion of AWE is more involved because the fully populated boundary integral 
submatrix of the system has a more complex dependence on frequency. In 
this case we may instead approximate the full submatrix with a spectral 
expansion of the exponential boundary integral kernel to facilitate the ex- 
traction of the system moments. This approach does increase the complexity 
in implementing AWE. However, AWE still remains far more efficient, in 
terms of CPU requirements when compared to the conventional approach of 
repeating the system solution at each frequency. 
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3 Numerical Implementation 

As an application of AW E to a full wave elect romagnet ir simulation, we con- 
sider t he evaluat ion of t he input impedance for a niicrost rip st uh shielded in a 
metallic rectangular cavity as shown in figure 1. The stub's input impedance 
is a strong function of frequency from 1 3 GHz and this example is therefore 
a good demonstration of AWE s capability. 

The shielded cavity is 2. .'58cm xfi.OOem x 1 ,06cm in size and t he niicrost rip 
stub resides on a 0.35m/ thick substrate having a dielectric constant of -5.2. 
The stub is 0.79m/ wide and A/2 long at 1.785 GHz and we note- that the 
back wall of the cavity is terminated by a metal-backed artificial absorber 
having relative constants of t T = (3.2. —3.2) and /i r = (1.0. — 1.0). 

As a reference solution, the frequency response of the shielded stub was 
first computed from 1 to 3 GHz at 40 MHz intervals (50 points) using a full 
wave finite element solution. To demonstrate the efficacy and accuracy of 
AWE we chose a single input impedance solution at 1.78 GHz in conjunct ion 
with the 4th order and 8th order AW E in (7) to approximate t he system 
response. Clearly the chosen number of poles or order of the expansion leads 
to different accuracies. As seen in Figure 2, the 4th order AWE representa- 
tion is in agreement with the real and reactive parts of the reference input, 
impedance solution over a 56% and 33% bandwidth, respectively. Surpris- 
ingly. the 8th order AWE representation recovers the reference solution over 
the entire 1-3 GHz band for both the real and reactive parts of the impedance. 
However, the CPU requirements for the 4th and 8th order approximations 
are nearly the same except for a few more matrix vector products needed for 
the higher order expansion. The number of these operations are of the same 
order as that of the AWE expansion and are much smaller than the size of 
the original numerical system. 

We conclude that the AWE representation is an extremely useful addition 
to electromagnetic simulation codes and packages for computing wideband 
frequency responses using only a few samples of the system solution. 


References 

[1] S. Kumashiro, R. Rohrer and A. Strojwas, “’Asymptotic waveform evalu- 
ation for transient analysis of 3-D interconnect structures,” IEEE Tran s. 


5 



( 'ompuh r-Aidrd Dfsign of Inf t prat id ("nruits and Sijstt ms. Yol. 12. 
No. 7. ])|> 988 990. Jul. 1993 


[2] L. Pillage and H. Rohrer, “AW F: Asymptotic Waveform Fst imat ion." 
7 ichn/cal Hr port. SRC-CMF Research (’enter for Computer-Aided De- 
sign. Carnegie-Mellon 1 Diversity. 1988 

[3] F. Chiprout and M. Xakhla. “Asymptotic* waveform evaluation and mo- 
ment matching for interconnect analysis. Norwell. kluwer Acad. Pubs. 
1 991 


[4] J. L. Volakis. A Chatterjee and J. Cong. " A class of hybrid finite ele- 
ment methods for electromagnetics: a. review,” fJIrctromagn. Warts 
Applications. Vol. 8, No. 9/10, pp. 1095-1124. 1994 

[5] J. Cong, J.L. Volakis. A.C. Woo and H. G. Wang, “A Hybrid Finite 
FI ement -Boundary Integral Method for the Analysis of Cavity-backed 
Antennas of Arbitrary Shape,’' IEEE Trans Antenna and Propagat .. 
Vol. 42, No. 9. pp 1233 1242, 1991 

[6] .Joseph Leliner, "Partial fraction decompositions and expansions of 
zero,” Trans. Amcr. Math. Sac.. Vol. 87, pp 130-143, 1958 


6 



45 



(a) 



Figure 3: Results of the 4th and 8th order AWE implementation using a 
single point expansion at 1.78 GHz. (a) Real part and (b) imaginary part of 
the Input impedance 
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